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The  combination  of  advanced  characterization  techniques  and  FEM-simulations  provided  detailed  infor¬ 
mation  about  losses  related  to  the  flowfield  geometry  of  a  metallic  interconnector  (MIC)  in  a  planar  SOFC 
repeat  unit.  The  presented  2-D  FEM  model  is  able  to  predict  the  repeat  unit  performance  decrease  due 
to  the  in  plane  ohmic  losses  in  the  electrodes  and  the  contact  resistance  between  electrode  and  MIC. 

The  performed  calculation  and  measurement  showed  an  increase  of  ohmic  losses  of  up  to  84%  when 
the  single  cell  was  contacted  with  a  MIC,  The  contact  resistance  adds  less  than  6%  on  the  cathode  side. 
The  in-plane  current  flow  from  the  contact  ribs  to  the  triple  phase  boundaries  under  the  gas  channels 
caused  94%  of  the  additional  ohmic  losses. 

Analysis  of  impedance  spectra  by  the  distribution  of  relaxation  times  and  a  subsequent  Complex  Non¬ 
linear  Least  Squares  fit  separated  gas  diffusion  from  the  total  polarization  losses.  Depending  on  flowfield 
design,  the  gas  diffusion  resistance  on  the  cathode  side  increased  up  to  +750%. 

For  high-performance  anode  supported  cells,  the  choice  of  cathode  flowfield  design  added  up  to  41% 
power  loss,  whereas  the  anodic  flowfield  design  was  of  minor  importance  (<1%  power  loss). 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  performance  of  high  temperature  solid  oxide  fuel  cells 
(SOFC)  at  the  single-cell  level  has  been  continuously  improved  over 
the  past  10  years  [1].  However,  at  the  stack  level,  the  power  den¬ 
sity  becomes  unprofitable  by  a  decrease  of  50%  [2].  A  search  for 
the  causes  leads  to  the  part  of  the  metallic  interconnector  (MIC). 
Apart  from  chromium  poisoning  of  the  cathode  [3-5],  additional 
losses  caused  by  the  MIC  geometry  [6]  reduce  the  cell  performance. 
Herein,  additional  ohmic  and  polarization  losses  are  created  by  the 
flowfield  geometry,  in  particular  (i)  non  ideal  contacting  through 
contact  resistance  or  limited  contact  area  and  (ii)  limited  gas  deliv¬ 
ery  through  gas  channels. 

Different  models  of  repeating  units  were  proposed  to  investi¬ 
gate  the  impact  of  MIC  design  on  the  performance  of  solid  oxide 
fuel  cells  (SOFC).  Nelson  et  al.  provided  a  continuum-level  model 
by  analyzing  the  effects  of  (i)  constriction  resistance,  (ii)  con¬ 
tacting  to  mass  transport,  and  (iii)  fuel  depletion  at  the  triple 
phase  boundaries  (TPB)  [7].  Tanner  and  Virkar  gave  an  analytical 
expression  for  the  ohmic  area-specific  resistance  (ASR0)  of  a  repeat 
unit.  They  used  differential  elements  considering  expressions  for 
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charge  transfer  resistance,  electrical  conductivities,  electrostatic 
potential,  and  MIC  geometry  [8].  A  time  dependent  model  was 
proposed  by  Gazzarri  and  Kessler  to  simulate  the  impedance 
spectra  of  a  repeat  unit.  They  developed  AC  and  DC  mass-  and 
charge  transport  equations  and  combined  them  with  the  gas  dif¬ 
fusion  and  the  concentration  polarization  within  the  electrodes 
[6]. 

Previously  published  modeling  attempts  [6-8]  present  complex 
physical  equations  with  a  large  quantity  of  parameters,  which  are 
difficult  to  determine  by  experiment.  In  contrast,  the  2D  FEM  model 
presented  here  is  based  on  experimental  data  extracted  from  (i) 
anode  supported  single  cell  performance  and  (ii)  electrical  con¬ 
ductivities  of  the  applied  materials.  The  2D  FEM  model  calculates 
(a)  the  current  density  and  (b)  the  electrical  potential  distribution 
of  a  planar  stack  repeat  unit  (MIC/contact/ASC/contact/MIC).  This 
numerical  model  can  be  adapted  to  other  cell  designs,  and  flowfield 
geometries. 


2.  FEM-model 

The  cross  section  of  applied  test  setup  shown  in  Fig.  1.  In  this 
study  we  investigated  ideal  contact  conditions  (a)  and  several  test 
set  ups  with  non  ideal  contact  conditions  of  a  repeat  unit  (b).  It  is 
obvious  that  test  setup  B  comes  closer  to  the  contact  conditions 
applied  in  a  stack.  Additional  losses  are  expected  here  due  to  (i)  the 
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Fig.  1.  Cross  sections  of  the  test  setups  with  an  ideal  contact  conditions  (a)  and  a  non  ideal  contact  conditions  (b).  The  single  repeat  unit  displayed  in  (b)  exhibited  a  contact 
rib  width  of  1  mm  and  a  gas  channel  width  of  4  mm. 


contact  resistance  between  electrode  and  flowfield,  (ii)  the  in  plane 
conduction  in  the  electrodes  as  well  as  (iii)  the  gas  transport  in  the 
thin  cathode  layer  beneath  the  contact  rib. 

It  has  to  be  emphasized,  that  the  transport  processes  (electronic 
and  ionic  current)  and  the  electrochemical  reactions  have  been  con¬ 
verted  into  a  purely  electric  steady  state  model.  Furthermore,  the 
in-plane  gas  diffusion  below  the  contact  ribs  is  not  considered.  The 
2D  model  was  implemented  in  the  finite  elemente  methode  soft¬ 
ware  COMSOL  Multiphysics  V3.5  (COMSOL  Inc.,  Burlington,  USA). 
The  discretisation  of  the  cell  and  contact  geometry  was  performed 
by  the  automatic  mesh  generator  of  COMSOL  without  local  mesh 
refining.  Altogether  5400  triangular  elements  were  generated.  The 
solving  process  of  a  current/voltage  (CV)  curve  (0. .  .2  A  cm-2,  step 
size  20  mA cm-2)  takes  not  more  than  45  s  in  a  conventional  work¬ 
station. 


2.1.  Geometry 

The  repeat  unit  shown  in  Fig.  1(b)  represents  the  geometric  con¬ 
tact  conditions  of  a  single  cell.  The  ASC  is  in  between  two  flowfields, 
comprises  of  both  gas  channels  and  contact  ribs  which  included  in 
the  interconnector  plates.  The  following  assumptions  are  made: 

(1)  the  contact  between  the  electrodes  and  the  contact  rib  are 
homogeneous  along  the  gas  flow  direction.  The  geometries  of  gas 
channels  run  parallel  to  the  contact  ribs  and  were  constant  in  the 
gas  flow  direction.  (2)  The  model  represents  an  infinite  cell  ele¬ 
ment  in  the  direction  of  the  gas  flow.  Therefore  the  gas  utilization 
can  be  neglected.  (3)  The  geometry  of  each  gas  channel  and  contact 
rib  is  equal.  As  a  consequence,  a  2D-model  considering  just  1/2  of 
the  contact  rib  (center  of  the  contact  rib  at  x  =  0)  and  1/2  of  the  gas 
channel  (center  of  the  gas  channel  at  x  =  L)  in  x-direction  describes 
the  impact  of  the  flowfield  design  on  the  cell  performance  (Fig.  2). 

Equations  and  boundary  conditions:  The  model  consists  of  4  sub 
domains:  the  anode  (anode  substrate  +  anode  functional  layer),  the 
two  layered  cathode,  and  one  MIC  per  electrode.  The  sub  domains 
are  implemented  as  electrical  conductive  sub-domains  (Fig.  2a). 
The  electric  charge  transport  within  the  domains  is  described  by 
the  2D  Laplace  equation  with  Neumann  boundary  condition  (no 
internal  species  generation): 

-V(<xd  Wd)  =  0  (1) 

where  ad  is  the  electronic  conductivity  and  Vd  the  electric  potential 
in  the  sub-domain.  The  boundary  conditions  are  defined  as  follows: 


(1)  The  current  is  impressed  as  a  current  source/sink  aty  =  0  and 
y  =  h.  The  electrical  potential  in  the  MIC  at  y  =  0  and  y  =  h  is  con¬ 
stant. 

(2)  Electric  isolation  for  the  symmetry  divisions,  i.e.  no  current  flow 

in  x-direction  at  x  =  0  and  x  =  L. 


Eq.  (2)  describe  the  total  current  flow  Jru  per  repeat  unit  and  Eq. 

(3)  is  used  to  calculate  the  total  voltage  loss  Uloss  per  repeat  unit. 


Im=[jydA  =  z[  Jcell(x)  \  dx 

(2) 

Ja  Jo  y 

fh 

Uloss  =  /  T(y)dy 

Jo 

(3) 

2.2.  Cell 

The  complex  electrochemical  reactions  and  the  transport  pro¬ 
cesses  of  an  ASC  are  represented  by  its  area  specific  resistance 
(ASR)  in  the  electrical  FEM-model.  The  open  circuit  voltage  (cor¬ 
responding  to  the  oxygen  activity  gradient  between  the  oxidant 
(air)  and  the  fuel  (H2,  5%  FI20))  is  not  considered  here  because  only 
the  voltage  losses  Uloss  of  the  repeat  unit  are  of  relevance.  Current 
voltage  ( CV)-characteristics  of  the  repeat  unit  ( UTU^ ),  are  simulated 
by  subtraction  of  Uloss  from  the  open  circuit  voltage  OCV  (Eq.  (4)): 

Uru(/)  =  OCV-Uloss(/)  (4) 

An  ASR-layer  is  inserted  in  between  the  electronically  conduct¬ 
ing  cathode  layer  and  the  electronically  conducting  anode  substrate 
(see  Fig.  2,  y  =  c).  The  ASR-layer  represents  the  ohmic  losses  of  the 
electrolyte  and  the  electrodes  as  well  as  the  polarization  losses 
of  the  two  electrodes  [9].  This  ASR-layer  is  implemented  as  a 
boundary  condition  in  the  COMSOL  FEM  model  describing  the  local 
current  density  jcell(x)|y=c  flowing  in  y-direction  through  the  ASR- 
layer  jcell(x)  as  a  function  of  the  local  electrical  potential  difference 
AU(x)|y=c  between  cathode  and  anode  at  y  =  c.  Both  jcell(x)|y=c  and 
AU(x)|y=c  vary  in  the  x-direction: 

AU(x)|y=c  =  Vcathode(x)|y=c  -  Vanode(x)|y=c  (5) 

To  parameterize  the  ASR-layer,  CV-characteristics  of  an  ideal 
contacted  cell  (see  Fig.  1 )  have  been  used.  A  homogeneous  current 
collection  and  gas  distribution  all  over  the  active  electrode  areas  is 
realized  by  a  gold  and  a  nickel  contact  mesh.  Under  these  ideal  con¬ 
ditions  all  transport  processes  (gas  and  bulk  diffusion,  electronic 
and  ionic  current  flow)  proceed  in  the  y-direction.  The  ohmic  resis¬ 
tance  of  the  anode  and  the  cathode  in  y-direction  was  calculated 
based  on  the  4-point  measured  conductivity  of  the  porous  electrode 
materials  and  the  thickness  of  the  layers.  Due  to  the  high  conductiv¬ 
ities  (anode:  150 103  Sm-1,  cathode:  current  collector  1300  Sm_1, 
functional  layer  320  S  m-1 )  and  the  low  thickness  (anode:  300  p,m, 
cathode:  75  p,m)  the  ohmic  resistance  is  fairly  below  0.2  m£2  cm1 2. 
In  the  end,  the  part  of  the  ohmic  losses  iny-direction  of  theses  layers 
in  comparison  with  the  sum  of  ohmic  resistance  of  the  electrolyte 
and  the  polarization  losses  of  the  electrodes  can  be  neglected  in 
case  of  the  ideal  contacting. 

The  total  amount  of  ohmic  and  polarization  losses  of  the  ideal 
contacted  cell  are  measured  by  CV-characteristics.  Considering  the 
ohmic  losses  in  the  electrodes  iny-direction  are  negligible,  the  CV- 
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(  14  rib  width  and  14  cannel  width): 


Fig.  2.  2D  FEM-Model  geometry:  the  repeat  unit  section  considered  in  the  model 
consists  of  1/2  of  the  contact  rib  and  1/2  of  the  gas  channel,  (a)  Implemented  sub 
domains,  (b)  mesh,  (c)  simulated  electrical  potential  distribution. 
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Fig.  3.  Test  geometry  and  calculated  surface  potential  of  the  cathode  for  differ¬ 
ent  ASRCc-values.  The  potential  difference  between  contact  and  surface  potential 
includes  the  contributions  of  the  contact  resistance  and  the  losses  caused  by  in  plane 
current  flow.  The  measured  values  match  with  an  ASRCC  of  3.4  mQ  cm2  (T  =  800  °C, 
P(02)cathode  =  0.21  atm,  p(H20)an0de  =  0.05  atm). 


characteristic  of  the  ideal  contacted  cell  provides  the  ASRcell  as  a 
function  of  the  electrical  potential  difference  All  at  y  =  c: 


ASRcell(At/)  = 


A  U 

Jcell(AU) 


OCV  —  UCell 
Jcell(OCV  -  Uce\\) 


(6) 


The  ASRcell  is  only  a  function  of  the  electrical  potential  difference 
AU  at  constant  temperature  and  gas  composition.  The  ASRcen(  A U) 
of  the  ideal  contacted  cell  was  implemented  in  the  boundary  con¬ 
dition  between  the  cathode  and  the  anode  domain  to  calculate  the 
local  current  density  in  y-direction  through  the  ASR-layer: 


Jcell(x)|y=c 


AuwU 

ASRceli(AU(x)|y=c) 


(7) 


2.3.  Contact  resistance 


The  contact  resistances  in  between  MIC  and  both  electrodes 
were  likewise  implemented  as  boundary  conditions  in  between 
the  concerned  domains.  In  this  study  the  area  specific  contact 
resistance  in  between  gold-flowfield  and  cathode  was  directly  eval¬ 
uated  from  the  repeat  unit  measurements.  The  simulated  electrical 
potential  distribution  in  the  cathode  was  fitted  to  the  electri¬ 
cal  potential  distribution  measured  by  means  of  potential  probes 
(Fig.  3)  using  the  cathode  contact  resistance  ASRCC  as  the  free  fit 
parameter. 

Total  ohmic  resistance:  To  calculate  the  total  ohmic  resistance 
of  the  repeat  unit  (ASR0(ru))  the  boundary  condition  (Eq.  (7))  in 
between  the  cathode  and  anode  sub  domains  were  modified.  The 
expression  ASRcen(  A U)  (Eq.  (6))  is  replaced  by  a  constant  resistance 
value  which  give  the  ohmic  resistance  of  the  electrolyte  (Eq.  (8)). 
The  ASR0  of  the  ideal  contacted  cell  (ASR0(ideai))  were  determined 
by  electrochemical  impedance  spectroscopy  (EIS)  and  provide  the 
ohmic  resistance  of  the  electrolyte  and  the  negligible  ohmic  con¬ 
tribution  of  the  electrodes  as  described  above.  The  calculation  of 
ASRo(ru)  is  carried  out  by  the  first  Ohm’s  law  (Eq.  (9)): 


Jcell(x)|y=c 


At/(*>U 

AS  Rq(  ideal) 


(8) 


ASRo(ru)  = 


U\oss 

7(ru) 


(9) 


3.  Experimental 


The  ASCs  within  this  study  are  based  on  50  mm  x  50  mm  anode 
supported  cells.  The  cosintered  porous  anode  substrate  (Ni/8YSZ- 
cermet),  the  anode  functional  layer  (Ni/8YSZ-cermet)  and  the  gas 
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Fig.  4.  Cross  sections  of  the  five  different  contact  geometries  applied  in  this  study  (a)  ideal  contact  (provides  ideal  performance  data),  (b)  test  geometry  (to  validate  the 
model)  and  (c),  (d),  (e)  flowfields  configuration  (to  investigate  the  limiting  factors  on  performance  of  single  cell). 


tight  electrolyte  (8YSZ)  were  manufactured  by  CeramTec  AG  (Mark- 
tredwitz)  [10].  Subsequently  a  double-layered  La0.65Sr0.3MnO3 
(LSM,  functional  layer)/LSM-8YSZ  (current  collector)  cathode  was 
applied  by  screen  printing  and  sintering.  Details  of  the  manufactur¬ 
ing  procedures  of  the  LSM-cathode  can  be  found  in  [12].  The  active 
area  of  the  working  cathode  was  1 0  mm  x  1 0  mm  and  two  auxiliary 
(OCV  probe)  electrodes  in  gas  flow  direction  in  front  of  and  behind 
the  cathode  [9].  ASCs  and  repeat  units  were  tested  in  a  single  cell 
test  bench  [13]  and  characterized  by  CV-  and  EIS  measurements. 

Five  different  contact  configurations  were  applied  in  this  study, 
as  displayed  in  Fig.  4:  (a)  ideal  contact  by  a  gold-mesh  (cathode)  and 
a  nickel-mesh  (anode)  providing  a  homogeneous  current  collection 
and  gas  distribution,  (b)  a  test  geometry  and  (c),  (d),  (e)  two  flow- 
field  designs  distinguishable  by  contact  rib  width  and  gas  channel 
width.  The  flowfields  were  made  of  gold  to  avoid  corrosion  effects 
at  the  cathode  and  nickel  for  the  anode  side,  respectively. 

Fig.  4(a)  shows  the  ideal  contacted  single  cell.  This  setup 
provides  the  ideal  performance  data  (ASRcell(AD))  which  was 
implemented  in  the  model.  The  model  was  validated  with  the  test 
geometry  in  Fig.  4b.  The  cathode  was  contacted  in  gas  flow  direc¬ 
tion  by  a  gold  wire  (1  mm  diameter,  10  mm  length)  at  x  =  5mm. 
The  contact  on  the  cathode  side  causes  an  in-plane  current  flow  in 
the  cathode  in  x-direction.  The  resulting  voltage  distribution  was 
measured  by  two  potential  probes  contacting  the  cathode  surface 
at  Xi  =4.2  mm  and  x2  =  7.8mm,  respectively  (see  also  Fig.  3).  The 
setups  in  Fig.  4(c)-(e)  were  used  to  investigate  the  cell  performance 
depending  on  the  flowfield  geometry:  (i)  the  anode  was  contacted 
by  a  nickel  flowfield  and  the  cathode  was  still  ideal  contacted  (gold- 
mesh).  The  flowfield  had  a  rib  width  of  1  mm  and  a  channel  width 
of  4  mm  (see  Fig.  4  c),  (ii)  the  cathode  was  contacted  by  a  gold  flow- 
field  and  the  anode  by  a  nickel  mesh.  The  dimension  of  the  flowfield 
is  rib  width  of  1  mm  and  a  channel  width  of  4  mm  (Fig.  4d).  (iii)  The 
cathode  was  contacted  by  a  gold  flowfield  and  the  anode  by  a  nickel 


mesh.  The  dimension  of  the  flowfield  is  a  channel  width  of  2.6  mm 
and  rib  width  of  2.7  mm,  (Fig.  4e).  The  electrical  potential  was 
measured  by  potential  probes,  a  spherical  (~1  mm  diameter)  gold 
(cathode),  respectively  platinum  (anode)  contact,  pressed  against 
the  electrode  surface  at  a  defined  x-position  (Figs.  4b-e  and  3).  The 
electrical  potential  at  the  electrode  surface  was  evaluated  by  mea¬ 
suring  the  potential  difference  vs.  a  reference  point,  i.e.  the  gold 
respectively  the  nickel  flowfield.  The  measured  potential  gives  the 
sum  of  the  voltage  losses  through  the  contact  resistance  and  the 
voltage  losses  causes  by  the  in  plane  current  flow  (x-direction). 

The  cells  were  operated  under  ambient  pressure  with  different 
N2/O2  mixtures  at  the  cathode  side  and  varying  H20/H2  mixtures 
at  the  anode  side.  The  total  anodic  and  cathodic  gas  flow  rates 
were  maintained  at  a  constant  value  of  0.25  ml  min-1  during  all 
experiments,  resulting  in  a  fuel  utilization  of  6%  and  an  oxidant 
utilization  of  12.6%  at  a  current  density  of  2  A  cm-2.  Additionally  a 
mixture  of  34%  H2 , 45%  N2 ,  and  2 1  %  H2  O  by  a  lower  gas  flow  rate  of 
199  ml  min-1  were  applied  to  investigate  more  realistic  operation 
condition  at  the  anode  side.  The  cells  were  tested  at  an  operat¬ 
ing  temperature  of  800  °C.  Potential  probes  measurements  and 
performance  data  were  carried  out  at  constant  current  load  and 
within  CV-characteristic  measurements.  Impedance  spectra  were 
measured  under  open-circuit  conditions  with  a  Solartron  1260  fre¬ 
quency  response  analyzer  in  a  frequency  range  from  100  mHz  to 
1  MHz.  The  amplitude  of  the  current  stimulus  was  adjusted  in  order 
to  achieve  a  voltage  response  <12  mV.  This  procedure  leads  to  a 
good  signal-to-noise  ratio  (SNR)  of  the  measured  data  and  fulfils 
the  condition  of  linearity  in  the  operating  point. 

4.  Results  and  discussion 

Impact  of  the  anode  flowfield:  The  electrical  potential  differ¬ 
ence  between  the  anode  substrate  and  the  anode  flowfield  was 


M.  Kornely  et  al.  /  Journal  of  Power  Sources  196  ( 201 1 )  7209-7216 


7213 


Fig.  5.  Measured  C/V-characteristics  of  applied  anode  contact  geometries  at  differ¬ 
ent  gas  composition:  (i)  H2  with  5%  H20,  and  (ii)  34%  H2,  45%  N2,  and  21%  H20 
(T=800°C). 

measured  for  two  different  anode  contact  geometries  applying 
an  ideal  cathode  contacting  (Fig.  4(a)  and  (c)).  The  voltage  loss 
related  to  the  anode  flowfield  was  always  below  1  mV.  The  ASR- 
contribution  of  the  in  plane  conduction  (anode  substrate  and  nickel 
flowfield)  and  the  contact  resistances  (anode  substrate/nickel- 
flowfield)  was  fairly  below  1  m^cm2.  This  is  in  good  agreement 
with  the  simulation  results  (Fig.  2c),  the  simulated  and  measured 
potential  difference  between  the  anode  substrate  and  the  anode 
flowfield  was  only  0.2  mV.  Due  to  its  high  electronic  conductivity 
the  electrical  potential  within  the  anode  substrate  and  the  flow- 
field  is  nearly  constant.  A  comparative  analysis  of  CV-characteristics 
(Fig.  5)  and  impedance  spectra  of  an  ideal  contacted  cell  (cathode: 
gold  mesh,  anode:  nickel  mesh)  and  of  a  cell  with  ideal  contacted 
cathode  (gold  mesh)  and  a  nickel  flow  field  at  the  anode  side 
showed  neither  any  difference  in  performance  nor  in  the  polariza¬ 
tion  (measured)  or  the  ohmic  (measured  and  calculated)  resistance 
(Table  1). 

These  measurements  were  carried  out  for  a  p(FI20)  in  the  fuel 
(H2)  of  5%  and  a  mixture  of  H2,  N2,  and  H20  to  consider  more  real¬ 
istic  stack  operation  conditions.  The  experimental  and  modelling 
results  showed  that  the  anode  flowfield  geometry  investigated  in 
this  study  has  a  negligible  impact  on  the  cell  performance  (Fig.  5). 

However,  an  anode  flowfield  with  a  larger  contact  rib  width 
would  further  decrease  the  ohmic  resistance  but  increase  the  gas 
diffusion  polarization.  Nevertheless,  this  would  not  improve  the 
performance  of  the  repeat  unit  significantly,  because  the  ohmic 
losses  are  already  close  to  zero  (<1^2 cm2).  These  calculations 
and  considerations  clearly  explain  why  no  further  experimental 
and  modelling  studies  were  carried  out  for  the  anode  flowfield 
design.  Fig.  5  shows  the  C/V-c haracteristics  of  ideal  contacted  cath¬ 
ode  and  both  anode  contact  designs  (ideal,  flowfield  design  1)  for 
a  H2/N2/H20  and  H2/H20  conditions,  respectively.  These  perfor¬ 
mance  measurements  confirm  that  the  impact  of  the  applied  anode 
flowfield  design  is  negligible  for  a  wide  range  of  fuel  gas  composi¬ 
tions.  Naturally,  the  anode  flowfield  design  has  to  be  reconsidered 
when  a  more  realistic  fuel  gas,  as  natural  gas  or  as  diesel  reformate 
will  be  examined. 

Application  of  potential  probes:  The  test  geometry  was  used  to 
measure  the  electrical  potential  distribution  at  the  cathode  surface. 
By  means  of  the  two  potential  probes  the  contact  resistance  ASRCC 
between  the  gold  contact  wire  and  the  cathode  was  identified.  The 
diagram  in  Fig.  3  shows  the  simulated  electrical  potential  Vcathode(x) 
of  the  cathode  surface  (y  =  c  +  75fxm)  for  different  ASRcc-values. 
Because  all  of  the  other  parameters  in  the  model  are  fixed,  the  con¬ 


Fig.  6.  Calculated  (lines)  and  measured  (symbols)  CV-characteristics  of  an  ideal 
contacted  cell  and  the  two  cathode  flowfield  designs  1  and  2  (T=800oC, 

p(02)cathode  =  0.21  atm,  p(H20)anode  =  0.05  atm). 

tact  resistance  ASRCC  can  be  evaluated  by  adapting  the  simulated 
electrical  potential  distribution  of  the  cathode  to  the  measured 
one.  In  case  of  a  contact  resistance  of  ASRCC  =  3.4  m£2  cm2  the  sim¬ 
ulated  electrical  potential  fully  agrees  with  the  measured  values  at 
Xi  =  4.2  mm  and  x2  =  7.8  mm.  The  application  of  a  potential  probe  in 
combination  with  the  FEM-model  enabled  us  to  evaluate  the  con¬ 
tact  resistance  between  electrode  and  flowfield  in  situ  during  the 
cell  measurements.  No  additional  contact  resistance  measurement 
is  required  to  determine  this  value. 

The  potential  difference  between  position  Xi  and  x2  at  the  cath¬ 
ode  surface  remains  unaltered  during  ASRCC  variation  (Fig.  3).  Thus, 
V2 1 ,  the  difference  in  electrical  potential  of  the  cathode  surface 
between  position  x !  and  x2,  is  independent  of  the  contact  resistance 
and  results  solely  from  a  voltage  drop  owing  to  the  in-plane  current 
flow  parallel  to  the  cathode/electrolyte  interface  (and,  therefore, 
also  between  Xi  and  x2).  This  in-plane  current  is  reduced  by  the 
current  flow  across  the  boundary  into  the  electrolyte.  Eq.  (7)  facili¬ 
tates  a  spatially  resolved  calculation  of  the  current  density  through 
the  electrolyte  as  a  function  of  the  boundary  condition  ASRcen  (Eq. 
(6)).  The  very  good  agreement  between  measured  and  simulated 
values  of  V2.\  shows  that  the  simulated  current  density  in  the  cath¬ 
ode  and,  thus,  the  model  as  a  whole  is  able  to  precisely  describe 
the  current  density  distribution  in  the  cathode  depending  on  the 
contact  geometry. 

Impact  of  the  cathode  flowfield  design:  In  Fig.  6  the  sim¬ 
ulated  and  the  measured  CV  characteristics  are  shown  for  an 
ideal  contacted  cell  as  well  as  for  the  two  analyzed  cathode 
flowfield  designs.  In  case  of  flowfield  design  1  the  measured 
performance  (Pdesigni(o.7V) =  1-03  Won-2)  is  decreased  by  21%,  as 
compared  to  the  ideal  contacted  cell  (Pideai(o.7V)  =  1.3  W cm-2).  The 
flowfield  design  2  caused  a  performance  decrease  of  up  to  41% 
(^design2(o.7V)  =  0.76Wcm-2).  Measured  and  simulated  values  for 
the  ideal  contacted  cell  are  in  excellent  agreement.  This  result  con¬ 
firms  that  the  model  correctly  simulates  the  CV  characteristics. 
Measured  and  simulated  values  for  the  flowfield  design  1  show 
a  rather  small  deviation  varying  with  current  density  (error  <  3%). 
At  low  current  densities  the  measured  values  are  slightly  above  the 
simulated  whereas  at  high  current  densities  the  simulated  voltage 
exceeds  the  measured.  Measured  and  simulated  values  for  the  flow- 
field  design  2  show  a  severe  deviation.  Moreover,  the  simulations 
result  in  a  superior  CV-characteristic  for  flowfield  design  2,  whereas 
the  measurements  show  a  superior  CV-characteristic  for  flowfield 
design  1.  The  deviation  between  the  model  and  the  experimental 
values  show  a  nonlinear  increase  with  increasing  current  density. 

Measured  and  calculated  ASR0(m)  are  shown  in  Table  1,  and  in 
Fig.  7.  The  values  are  in  good  agreement  for  the  ideal  contacted 
cell  and  both  cathode  flowfield  designs.  In  case  of  design  1  and  2 
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Table  1 

Calculated  and  measured  ohmic  resistance  ASRO  and  the  measured  polarization  resistance  ASRpol  of  an  ideal  contacted  cell,  the  anode  flowfield  design  1,  and  the  two  cathode 
flowfield  designs  1  and  2  at  OCV-condition  ( T  =  800 °C,  p(02)cathode  =  0.21  atm,  p(H20)anode  =  0.05  atm). 


Cathode 


Ideal  contact 
Anode  design  1 
Cathode  design  1 
Cathode  design  2 


Anode 

ASR0  [m£2cm2] 

ASRp0i  [m^cm2] 

Measured 

Calculated 

Measured 

34 

34 

316 

34 

34 

315 

63 

64 

318 

42 

42 

355 

the  contact  resistance  ASRCC  is  below  2  m£2  cm2.  The  total  increase 
of  ASR0  (design  1  =  87%,  flowfield  design  2  =  24%)  is  mainly  caused 
by  the  in-plane  current  flow  from  the  contact  ribs  to  the  triple 
phase  boundaries  (TPB)  under  the  gas  channels  (design  1=94%, 
design  2  =  87%).  Due  to  the  channel  width  the  ASR0(design  2)  is  lower 
than  ASR0(design  i)  whereas  the  performance  displayed  an  oppos¬ 
ing  dependency.  This  is  a  strong  evidence  for  a  further  performance 
limiting  factor  (Fig.  7). 

Fig.  8  shows  the  simulated  cathode  potential  distribution 
^cathode (*)  of  the  cathode  surface  (y  =  c  +  75  fxm)  and  the  related  cur¬ 
rent  density  distribution  jcell(X)  |y=c  for  different  current  loads  in  case 
of  flowfield  design  2.  Using  the  procedure  described  above,  an  ASRCC 
of  1  m^cm2  was  determined  for  a  current  density  of  20  mA cm-2. 
Because  the  ASRCC  does  not  depend  on  the  current  density,  this 
value  should  be  applicable  for  all  simulations. 

With  increasing  current  density  the  measured  value  of 
VCathode(x  =  5  mm)  slightly  exceeds  the  simulated  values.  This  indi¬ 
cates  that  the  simulated  current  density  in  the  gas  channel  is  too 
low  and/or  the  simulated  current  density  under  the  contact  ribs 
is  too  high.  Considering  that  a  constant  current  density  under 
the  contact  ribs  requires  a  homogeneous  supply  of  02  the  in¬ 
plane  gas  diffusion  under  the  contact  rib  has  to  be  considered. 
The  in-plane  gas  diffusion  path  to  the  TPB  under  the  contact  area 
(max.  1200  |xm)  is  significantly  longer  than  in  the  case  of  the  ideal 
contacted  cell  (only  vertical  gas  diffusion,  ~75  p,m)  and  causes 
additional  losses.  Therefore  the  assumption  of  negligible  in-plane 
gas  diffusion  is  a  probable  reason  for  the  deviations  described 
above. 

To  further  analyze  this,  impedance  spectra  at  different  oxi¬ 
dant  gas  compositions  (0.5-21%  02  in  N2)  were  measured.  The 
impedance  measurement  and  data  analysis  was  carried  out  accord¬ 
ing  to  the  method  described  in  of  our  previous  papers  [9,10].  As 
it  can  be  seen  in  the  distribution  of  relaxation  times  (DRT)  [11] 
(Fig.  9a))  the  gas  diffusion  process  PiC  and  the  oxygen  reduction 
process  P2c  increase  with  decreasing  p(02)cath0de-  The  gas  diffu¬ 
sion  resistance  RD(cathode)  was  obtained  for  low  p(02)cathode  (<2%) 
by  a  CNLS  fit  using  the  equivalent  circuit  depicted  in  Fig.  9b) 
[9].  A  direct  evaluation  of  RD(cathode)  fr°m  an  impedance  spectra 


measured  in  air/H2  (5%  H20)  is  impossible  because  the  small  gas 
diffusion  process  PiC  at  the  cathode  (RD(cathode)(air)<  5  m^2  cm2, 
relaxation  frequency  ~9Hz)  is  hidden  by  the  gas  diffusion  process 
P1A  in  the  anode  substrate  (R1A  =  130  m£2  cm2,  relaxation  frequency 
~4  Hz).  By  well  selected  operating  conditions,  namely  a  high  p(H20) 
at  the  anode  (63%  H20)  and  a  low  p(02)  at  the  cathode  (0.4 . . .  8%), 
RiA  is  decreased  to  30  m^  cm2,  Rotcathode)  increases  and  the  relax¬ 
ation  frequency  of  PiC  is  shifted  to  lower  values.  In  this  case  the 
resistance  RiA  was  kept  fixed  during  the  entire  fit  procedures.  Eq. 
(10)  was  used  to  calculate  the  RD(cathode)  in  air  [14].  This  equation 
describe  the  gas  diffusion  resistance  in  porous  electrodes  for  a  per¬ 
fect  one-dimensional  diffusion  (y-direction)  in  dependency  of  the 
p(02)cathode>  th e  temperature,  and  structure  parameters. 

In  Fig.  9(c)  Ro(cathode)  is  displayed  for  the  ideal  contacted  cell 
and  the  two  flowfield  designs.  It  is  obvious  that  the  RD(cathode)- 
values  and  their  p(02 ^dependencies  are  significantly  influenced 
by  the  flowfield  design.  The  gas  diffusion  resistance  increases 
with  increasing  contact  rib  width.  Due  to  the  increase  in  plane 
gas  diffusion  under  the  wide  contact  ribs  (x-direction)  the  p(02)- 
dependence  could  not  further  describe  by  Eq.  (10)  [14],  which 
assumes  a  one-dimensional  gas  diffusion.  Due  to  this  the  Ro(cathode) 
in  air  in  case  of  wide  ribs  was  calculated  by  a  linear  extrapolation 
of  the  fitted  values. 

-  (jj)  (p(02)clt„  _1) 

(10,33*,o5S5;)  '  (,0> 

The  change  in  p(02)-dependence,  which  was  observed  in  case 
of  the  flowfield  design  2  (wide  contact  ribs),  has  to  be  attributed 
to  a  coupling  of  the  in  plane  gas  diffusion  resistance  in  the  cathode 
layer  and  the  area  specific  resistance  of  the  cell. 

These  results  indicate  for  the  cathode  side,  that  the  in  plane 
gas  diffusion  below  the  contact  rib  must  be  considered.  Obviously, 
the  FEM-model  of  the  repeat  unit,  which  does  not  include  in-plane 
gas  diffusion  under  the  contact  ribs,  is  insufficient  for  a  contact  rib 


Fig.  7.  Calculated  ASR0(rU)  (filled  symbols)  and  measured  EIS  spectra  of  an  ideal  contacted  cell  and  the  two  cathode  flowfield  designs  1  and  2  at  OCV-condition  (T=  800  °C, 
P(02)Cathode  =  0.21  atm,  p(H2O)anOde  =  0.05  atm). 
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Fig.  8.  Simulated  potential  Vcathode(x)  (a)  and  current  density  jceu(x)\y=c  (b)  distribution  for  different  currents  Jru  ranging  from  100  mA  to  0.8  A  for  cathode  flowfield  design  2. 
The  measured  values  Vcathode  (x  =  5mm)  are  displayed  as  filled  symbols  (T=800°C,  p(02)cathode  =  0.21  atm,  p(H20)anode  =  0.05  atm). 
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Fig.  9.  (a)  Distribution  of  relaxation  times  corresponding  to  measured  impedance  data  (T=  800  °C,  OCV-condition,  p02(cathode)  =  0.01 1 , . . .,  0.21  atm,  pH20(anode)  =  0.05  atm). 
For  cathode  flowfield  design  1,  (b)  equivalent  circuit  model  with  cathode  gas  diffusion  process  P1C  highlighted,  and  c)  gas  diffusion  resistance  RD(cathode)(p(02)cathode)  for  the 
ideal  contacted  cell  and  the  two  cathode  flowfield  designs. 


width  exceeding  1  mm.  In  this  case,  in-plane  gas  diffusion  in  the 
cathode  layer  must  be  integrated  in  the  FEM  model. 

In  Fig.  10  the  ohmic  resistance  ASR0  and  the  gas  diffusion  resis¬ 
tance  Ro(cathode)  is  displayed  for  the  ideal  contacted  and  the  two 
flowfield  designs.  The  increase  in  ASR0  and  Ro(cathode)*  which  is 
related  to  the  electrical  conduction  beneath  the  gas  channel  and 
the  gas  diffusion  beneath  the  contact  ribs,  respectively,  correspond 
with  the  different  contact  rib  and  gas  channel  width  of  the  2  cathode 
flowfield  designs. 

Our  experiments  and  the  FEM-model  show  an  increase  of  ohmic 
losses  of  84%,  when  the  cathode  is  contacted  by  a  MIC  with  narrow 
contact  ribs  (design  1 :  1  mm  contact  rib  width,  4  mm  gas  channel 
width).  The  increase  in  ohmic  losses  is  narrowed  down  to  +24%, 
when  the  cathode  is  contacted  by  a  MIC  with  wide  contact  ribs 
(design  2:  2.4  mm  contact  rib  width,  2.6  mm  gas  channel  width). 
This  observation  cannot  be  explained  by  the  contact  resistance 
between  MIC  and  cathode,  as  it  was  less  than  3%  (3.4 m£2 cm2). 
Fact  is  that  the  in-plane  current  flow  from  the  contact  ribs  to  the 
triple  phase  boundaries  is  responsible  for  these  ohmic  contribu¬ 
tions.  As  the  MIC  with  wide  contact  ribs  (design  2)  has  a  smaller 
gas  channel  width  (2.6  mm),  the  in-plane  ohmic  losses  in  the  cath¬ 
ode  are  smaller.  This  part  of  the  ohmic  losses  is  governed  by  (i) 


Fig.  10.  Increase  of  the  ohmic  resistance  ASR0  and  the  gas  diffusion  resistance 
Rd( cathode)  for  the  two  flowfield  designs.  In  case  of  design  2  the  RD(cathode)  increase 
about  750%  (T  =  800  °C,  OCV  condition  p02(cathode)  =  0.21  atm,  pH20(anode)  =  0.63 
atm). 
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thickness  and  (ii)  electronic  conductivity  of  the  cathode.  As  a  con¬ 
sequence,  the  performance  of  the  repeat  unit  can  be  increased  by 
layout  and/or  materials  choice. 

Based  on  our  experiments  at  OCV  conditions,  the  gas  diffusion 
resistance  increased  by  48%  for  the  flowfield  design  with  narrow 
contact  ribs  (design  1:  1  mm  contact  rib  width,  4  mm  gas  channel 
width).  As  the  ideal  contacted  cell  exhibits  a  diffusion  resistance  of 
4.4m^cm2,  this  increase  would  correspond  to  6.5  m^  cm2  only. 
In  contrast  to  this,  in  case  of  wide  contact  ribs  (design  2:  2.4  mm 
contact  rib  width,  2.6  mm  gas  channel  width)  the  increased  dif¬ 
fusion  resistance  (+750%)  leads  to  a  remarkably  high  value  of 
37.5m£2cm2. 

5.  Conclusions 

The  combination  of  advanced  characterization  techniques 
and  FEM-simulations  provided  detailed  information  about  losses 
related  to  the  flowfield  geometry  of  a  MIC  in  a  planar  SOFC  repeat 
unit.  The  2D  FEM  model  developed  in  this  study  is  able  to  predict 
the  repeat  unit  performance  decrease  due  to  the  in  plane  ohmic 
losses  in  the  electrodes  and  the  contact  resistance  between  elec¬ 
trode  and  interconnector.  Appropriate  analysis  of  impedance  data 
in  combined  with  a  calculation  of  the  distribution  of  relaxation 
times  (DRT)  separated  the  gas  diffusion  losses  below  the  contact 
ribs  from  the  total  polarization  losses.  Experimental  and  modeling 
results  showed  that  the  flowfield  design  on  the  cathode  side  had  a 
large  effect  on  the  repeat  unit  performance  whereas  the  flowfield 
design  on  the  anode  did  not  influence  the  performance  under  the 
testing  conditions  applied  in  this  study. 

The  higher  resistance  of  a  repeat  unit  has  to  be  attributed  to  in¬ 
plane  transport  processes  in  the  cathode  layer:  (i)  in-plane  ohmic 
losses  underneath  the  gas  channels,  which  increase  with  gas  chan¬ 
nel  width  and  (ii)  in-plane  gas  diffusion  losses  underneath  the 
contact  ribs  which  increase  with  contact  rib  width.  As  a  conse¬ 
quence,  the  best  flowfield  design  was  identified  as  a  combination 
of  a  narrow  gas  channel  width  a  narrow  contact  rib  width. 


With  respect  to  high-performance  anode  supported  cells,  the 
flowfield  design  on  the  cathode  side  adds  up  to  41%  power  loss, 
whereas  the  flowfield  design  on  the  anode  side  is  of  minor  impor¬ 
tance  (<1%  power  loss). 
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